UNCLASSIFIED 


END 

DATE 


9-77 


WISCONSIN  UNIV  MADISON  MATHEMATICS  RESEARCH  CENTER 
nonneutral  plasma  expansion. (U) 

MAY  77  K E LONNGREN  DAAG29-75-C-0024 


F/G  20/9 


MRC-TSR-1747 


Mathematics  Research  Center 
University  of  Wisconsin— Madison 
610  Walnut  Street 
Madison,  Wisconsin  53706 


O 


May  1977 


(Received  April  25,  1977) 


Sponsored  by 

U.S.  Army  Research  Office 
P.O.  Box  12211 
Research  Triangle  Park 
orth  Carolina  27709 


Approved  for  public  release 
Distribution  unlimited 


National  Science  Foundation 
Washington,  D.C.  20550 


UNIVERSITY  or  WISCONSIN  - MADISON 
MATHEMATICS  RESEARCH  CENTER 


NONNEUTRAL  PLASMA  EXPANSION 

* 

K,  E.  Lonugren 


Technical  Summary  Report  #1747 
May  1977 

ABSTRACT 

The  self  similar  expansion  of  nonneutral  plasmas  is 
examined . using  various  models.  Analytical  solutions  are 
obtained  for:  (1)  The  expansion  of  a Maxwellian  electron 

cloud  governed  by  Ohm's  Law;  (2)  The  expansion  of  an  electron 
cloud  using  a mobility  model;  (3)  The  expansion  of  an  electron 
cloud  using  a cold  plasma  model ^ and  (4)  The  expansion  of  a 
charge  particle  cloud  with  a temporally  decaying  nonlinear 
diffusion  coefficient. 
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NONNEUTRAL  PLASMA  EXPANSION 

* 

K.  E,  lonngren 

In  a previous  report/^vre  presented  a methodical  procedure  for  obtaining  the  self 
similar  variables  using  some  of  the  techniques  of  Lie  Group  Theory.  General  physical 
phenomena  which  occur  in  plasmas  and  could  be  modeled  with  equations  amenable  to  a self- 
similar treatment  were  presented  along  with  their  solutions.  An  extensive  list  of  ref- 
erences was  also  given.  The  purpose  of  this  report  is  to  summarize  a pot-pourri  of 
further  examples  which  fall  under  the  umbrage  of  "Nonneutral  Plasma  Expansion". 

With  the  increased  interest  in  relativistic  electron  beam  fusion  devices,  consider- 
able attention  has  been  given  to  the  study  of  nonneutral  plasmas.  The  recent  monograph 
(2) 

by  Davidson  summarizes  the  past  work  and  motivates  our  interest  in  examining  various 
aspects  of  this  problem.  In  this  report,  we  shall  examine  three  sets  of  fluid  equation 
models  and  using  the  technique  of  "Self  Similar  Solution  of  Partial  Differential  Equations 
obtain  analytical  solutions. 

In  addition,  a recent  experiment  on  the  Wisconsin  Multipole  suggests  that  the  cross- 
field diffusion  coefficient  can  be  modeled  with  a one-dimensional  diffusion  equation 
where  the  diffusion  coefficient  is  nonlinear  and  exponentially  decaying  in  time. A 
self-similar  solution  of  this  problem  shall  also  be  given. 

In  Section  II,  we  describe  the  various  models  and  present  the  self-similar  solutions 
The  models  are:  (1)  The  expansion  of  a Maxwellian  electron  cloud  governed  by  Ohm's  Law; 

(2)  The  expansion  of  an  electron  cloud  using  a mobility  model;  (3)  The  expansion  of  an 
electron  cloud  using  a cold  plasma  model;  and  (4)  The  expansion  of  a charge  particle 
cloud  with  a temporally  decaying  nonlinear  diffusion  coefficient.  Section  III  is  the 
conclusion. 


On  leave  from  the  University  of  Iowa. 
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II . Various  Problems 

To  analyze  the  problems,  we  follow  the  procedure  given  in  reference  1 and  only 
describe  the  physical  phenomena,  list  the  PDE,  the  self  similar  variables,  the  ODE  and 
the  solution  to  the  ODE  without  repeating  the  details  of  the  procedure  in  each  case. 

(4) 

1)  Expansion  of  a Maxwellian  electron  cloud  governed  by  Ohm's  Law 

In  this  example,  we  examine  the  expansion  into  a vacuum  of  a thermalized  electron 
cloud,  described  by  an  isothermal  Maxwellian  distribution.  It  is  assumed  that  collisions 
are  sufficiently  frequent  such  that  we  can  speak  of  a conductivity  for  the  medium. 

The  governing  equations  are: 

I)  Equation  of  continuity 

P,  = 0 (1) 

II)  Ohm's  Law 

j = and  (2) 


III)  Maxwellian  electrons 

P = n^q  expIq'J’/kgT^l 


(3) 


where  all  symbols  are  standard.  By  differentiating  (3)  with  respect  to  x and  substitut- 
ing (2) , we  obtain 


j 


k„T 
B e 


1 

P 


(4) 


In  normalized  units,  (4)  and  (1)  lead  to 


where 


(y  •I'  i - = 0 

<ii  y y T 


(5) 


'I'  = P/n^q,  y = ^ “ t/(e^/o) 

and  is  the  Debye  length  and  ® relaxation  time. 

We  find  that  (5)  admits  a self-similar  solution  of  the  form 

ri(t)  = Tl/)(y,T),  £;  = y/T  . (6) 

The  solution  satisfies  the  conservation  law  that  electron  charge  is  conserved  in  space. 
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The  boundary  conditions  which  are  germane  to  this  problem  are: 


ll. 


a)  iHo.t)  - i|<u/T  ■*  n(C=o)  = 


b)  iK“,T)  ■■  0 = 4'(y>0)  • = 0 i.e.  "consolidation". 


Substituting  (6)  into  (5),  we  obtain 


The  first  integral  of  (8)  is 


-Cn  * 


We  shall  further  impose  the  condition  that  the  current  j 0 as  y 


The  integral  of  (9)  is 


r/2  + c. 


From  (4) , 


this  transforms  to  — ri^  0 as  ^ . Using  this  and  (7b)  , Cj^  = 0 . 


The  constant  is  determined  from  (7a)  to  be  = I/’J'q  • terms  of  y and  T , the 


final  result  is 


. P 1 

<|<  = — — = 

V T(yV2T^  + l/lJ)jj) 


Charge  is  conserved  as  is  shown  below. 


/ ay  = / 

^ n n ■'  n n J 


0 "o'*  "o^  0 lC^/2  + l/i|;^) 


— = const. 


In  conclusion,  we  have  examined  the  expansion  of  an  electron  cloud  in  a vacuum. 


Under  conditions  where  Coulomb  forces  can  be  neglected,  this  calculation  could  model  an 


electron  cloud  expansion  in  a plasma  in  time  scales  short  with  respect  to  ion  motion. 


2)  Expansion  of  an  electron  cloud  using  a mobility  model 


Recently,  considerable  attention  has  been  given  to  the  problem  of  tlie  transient  be- 


havior of  the  bulk  electric  field  cUid  space  charge  distribution  in  semi-conductors  and  in 


the  conduction  in  dielectric  and  insulating  fluids.  It  has  been  found  prudent  to  use  a 


J 


mobility  model  where  the  velocities  of  charge  carriers  injected  from  an  emitting  electrode 
are  proportional  to  the  electric  field  through  their  mobilities  and  the  electric  field  is 
related  to  the  charge  densities  of  the  carriers  through  Gauss's  law.  Many  and  Rakavy 
and  Helfrich  and  Nark  were  probably  the  first  to  suggest  that  the  problem  could  be 
modeled  with  the  set  of  dimensionless  equations 


E = p 

X 


i + p.  = 0 
X t 


i = pE 


which  arc  Poisson's  equation,  the  equation  of  continuity  and  a mobility  definition  for 
current  respectively.  The  subscript  x and  t denote  a partial  differentiation  with 


respect  to  space  and  time. 

In  their  original  paper.  Many  and  Rakavy  obtained  a solution  to  (13)  by  looking 
for  the  characteristics  of  the  problem.  Subsequently,  this  approach  was  extended  by  Batra, 
Schechtman  and  Seki,^^^  Zahn,  Tsang  and  Pao, de  Oliveira  and  lerreira and  others. 

An  extensive  list  of  relevant  exi>erimental  observations  is  given  in  reference  9. 

As  the  problem  is  extremely  important,  we  suggest  an  alternative  technique  for  solu- 
tion which  will  describe  the  spatial  and  temporal  evolution  of:  I)  a fixed  electric  field 
and  II)  a constant  source  of  current  which  are  both  governed  by  (13) . The  technique  that 
we  shall  apply  is  to  find  the  "self-similar  solution"  of  this  set  of  partial  differential 
equations. 

Equation  (13)  can  be  written  as 

EE  + (E  )^  4 E ^ = 0 . (14) 

XX  X xt 


The  self  similar  variables  are 


— and  f,  = 


where  a/y  ajid  (1/y  arc  constants  wliich  will  be  specified  by  invariance  and  conserva- 
tion requirements.  Invariance  specifics  tlial  a - (<  = -y  . 


Substituting  (IS)  in  (14),  we  write 
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We  shall  obtain  solutions  for  (16)  subject  to;  I)  A fixed  electric  field  and  IT.;  a con- 
stant current  at  x = 0 requirement. 

I)  Electric  field  is  constant  at  x = 0 . We  choose: 


B = Y . 

With  these  values,  C = x/t  and  4*  = E . The  integral  of  (16)  with  these  constants  is 

(1-3 

tdiere  k^  is  a constant  of  integration.  The  constant  is  set  equal  to  zero  since 
E(x  = 0,t)  = 0 in  order  to  satisfy  space  charge  limited  conditions.  This  specifies 
♦(C  « 0)  =0  . The  solution  of  (17)  is  9 = C from  which  we  compute  that 


where  the  constants  x^  and  t^  have  been  explicitely  included  since  (13)  is  invariant 
to  translation. 

II)  Current  is  constant  at  x = 0 . We  choose: 

2a  ^ $ 

y " y ' 

with  this  choice,  we  have  5 a^id  9 “ — where  9 satisfies  (16)  which  becomes 

(9  - 209^^  - 9^  + (9^)^  = 0 . (20 


This  can  be  integrated  once  to 


(4.  - 20^^  + 4'  = . 


(21) 


The  constant  of  integration  is  set  equal  to  zero  since  we  require  that  E(x  = 0,t) 

be  zero  for  space  charge  limited  conditions  which  specifies  = 0)  =0  . The  integral 

of  (21)  is 

♦ = C (22) 


from  which  we  compute  that 


E = 


x+x 


0 


t+t 


0 


P 





x+x„ 


(t+to)' 


(23) 


where  again  tlie  constants  x^  and  t^  have  been  introduced  since  (13)  is  invariant  to 
translation.  Note  that  this  is  identical  to  (18) . 

In  (21) , we  can  also  obtain  the  solution  for  a non-space  charge  limited  condition 
(E(x  = 0,t)  ^ 0)  by  setting  the  constant  of  integration  equal  to,  say,  2i^  . The 

integral  of  (21)  can  then  be  written  as 

^ = C + <24) 

from  which  we  compute  that 


E 


x+x. 


(t  + t^)(i^  + 


(t+tp)' 


P = 


1 


(25) 


x+x. 


(t+t^) 


2 ^0 


where  again  the  constants  x^  and  t^  have  been  reintroduced.  Note  that  (25)  reduces 
to  the  space  charge  limited  case  for  i^  = 0 . 

In  conclusion,  we  have  shown  tliat  tlie  set  of  equations  w)iich  describe  the  Transient 
Space  Charge  Limited  Current  Problem  admit  self-similar  solutions  for  two  physically  in- 
teresting boundary  conditions.  These  solutions  are  valid  in  the  initial  stages  befoi e 
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the  particles  reach  a second  electrode  which  may  be  placed  at  x *=  L . 

3,  Expansion  of  an  electron  cloud  using  a cold  plasma  model 

A model  which  can  describe  the  behavior  of  an  electron  cloud  expansion  into  a fixed 
ion  background  is  to  assume  that  Poisson's  equation  is  an  initial  condition.  The  electro- 
static approximation  for  the  V X b Maxwell  equation  and  the  continuity  equation  assure 
that  Poisson's  equation  is  satisfied  for  all  time. 

The  basic  equations  are: 

+ Pj,  = 0 

mv^  + mw^  = -eE  (26) 

V X B ~ 0 = - pv 


which  are  the  equations  of  continuity  and  motion  and  Maxwell's  equation  respectively. 

OO 

The  self  similar  variables  which  satisfy  the  conservation  law  that  / ndx  = constant 

0 


are 


2 2 

e = E;  N = nt  j u = v/t;  C = x/t  . (27) 

These  are  the  same  self  similar  variables  that  were  obtained  in  an  earlier  study  of  the  set 

(26)  where  the  Ansatz  that  Maxwell's  equation  could  replace  Poisson's  equation  had  not 
(12) 

been  made.  In  the  earlier  study,  it  was  not  a pedestrian  tas)c  to  integrate  the  ODE. 

Substituting  (27)  into  (26),  we  now  obtain  the  ODE: 


-2N  - 2^N^  + (NU)^  = 0 

U - 2CU^  + UU^  = -e  (28) 

2Ce^  + NU  = 0 . 

A solution  of  this  set  is 

U = 2^  , N = 2 and  E = -2C  (29) 

from  which  we  write  the  solution  of  (26)  using  (27)  and  (29)  as 


E ■=  - 


2x 
. 2 


(30) 


In  conclusion,  we  find  that  in  the  final  self  similar  solution,  the  density  is 


-7- 


independent  of  position  at  the  end  rather  than  making  it  an  a priori  assumption  in  the 
calculation  as  did  Gintsburg  who  treated  a similar  problem. 

4.  Expansion  of  charged  particles  \s-ith  a temporally  decaying  nonlinear  diffusion 
coefficient 

In  recent  experiments  on  the  Wisconsin  Multipole,  it  was  confirmed  that  the  cross- 

(3) 

field  diffusion  coefficient  depended  on  time  and  amplitude  as 

-«t 

D ~ -E . (31) 

Incorporating  this  in  the  one  dimensional  diffusion  equation,  we  obtain 


(32) 


where  all  constants  except  a have  been  suitaiily  normalized  away.  A change  of  variables 


-at, 
C J 


(33) 


transforms  (32)  to 


n = 


T 


(34) 


The  self  similar  treatment  of  (34)  is  straightforward,  at  least  for  tlie  case  where 

OO 

the  conservation  law  / ndx  - constant  is  valid.  Tlie  self  similar  variables  are 

0 


V3  j r / 2/3 
N = ni  and  ^ = x/t 


(35) 


and  the  resulting  ODE  is 


-f(CN)^ 


[A 


(36) 


If  the  burst  of  particles  is  siTimetric  at  x = 0 such  that  n 


= 0,  and 


n(x=0,T)  = 1 


-2/3 


'x=0 


(36)  can  bo  integratc.d  twice  to  yield 


/~N  = 


2 + C /3 


(37) 


Using  (33)  and  (35)  in  (37),  we  finally  obtain 


n(x,t)  = 


4(i(l-e‘”''))2 


I2(i(l-e'«S,V3 


(38) 


We  note  that  (38)  gives  a reasonably  accurate  qualitative  description  to  the  ex- 
perimental results.  This  seems  true  even  though  the  time  scale  in  the  experiment  is 
long  such  that  normal  modes  have  been  exited. 

III.  Conclusion 

The  self  similar  behavior  of  four  plasma  phenomena  have  been  described. 
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